﻿ W şi Rao O INTRODUCERE ÎN DETECȚIA SEMNALELOR Problemele de detecție apar în: RADAR, telecomunicații, sonar, medicină, prelucrarea și recunoașterea imaginilor, seimologie automatizări ș a În toate aceste domenii trebuie să decidem dacă un eveniment este prezent sau nu și apoi să obținem mai multă informație despre el Detectarea evenimentelor ține de “teoria detecției”, cunoscută și sub denumirea de “testarea ipotezelor” sau “teoria deciziilor statistice” În anumite cazuri trebuie să decidem între două ipoteze posibile, situație în care se vorbește despre ipoteze binare Acesta este cazul detecției unei ținte prin RADAR Alteori trebuie să decidem între mai multe ipoteze, situație în care se vorbește despre ipoteze multiple Astfel de cazuri apar în telecomunicații, în clasificarea și recunoașterea formelor (pattern recognition) ș a Trebuie să decidem între două sau mai multe ipoteze, având la dispoziție date observate, care, prin natura lor sunt afectate de zgomot Datele observate fiind aleatoare, pentru detecție este necesară o abordare statistică Admitem că, uneori, deciziile luate pot fi eronate Urmărim ca probabilitatea deciziilor corecte să fie cât mai mare, iar rata erorilor cât mai mică sunt posibile două clase de detecție, cea bazată pe abordarea Neyman- Vom vedea că Pearson și cea bazată pe abordarea Bayes 1 Pentru exempificare, considerăm problema recepției unui semnal util, s(t), din care se prelevează un eșantion, x Ca urmare a prezenței unui zgomot aditiv de putere cunoscută eșantionul poate fi 0: x =0+w ; este prezent doar zgomotulH 1: x =A+w =1+w ; este prezent semnalul util de 1Volt, afectat de zgomotH 2 w 0,σ ∼N () Avem acces experimental doar la eșantionul x Cunoscând doar acest eșantion și caracteristicile statistice ale eșantionului de zgomot alb, gaussian, trebuie să spunem care dintre cele două ipoteze este adevărată Se poate vorbi și despre detectarea nivelului continuu de A=1 Volt din zgomot Deoarece zgomotul este de medie nulă, “bunul simț tehnic” ne spune că e bine să comparăm eșantionul x cu un prag A/2=0 5 Volt Decidem că: rată≥H Volt, semnalul A=1Volt este prezent sau 0x 0 5 este ipoteza adevă ratăH 0 5 Volt, semnalul A=1Volt nu este prezent sau 1x ⇒HHH () () 010Dacă Notăm PHHH; probabilitatea de a decide că ipoteza este adevărată, {} jii atunci când, de fapt, ipoteza este cea adevărată jH și rezultă că putem avea două tipuri de erori = xx () () pxH 0; Pragul testului, γ, este soluția ecuației γ===>xx xxH ()() {} FAPp d I Lα 0; ; unde : ∫ I *** L(x) se numește raport de plauzibilitae și este raportul a două plauzibilități, deoarece se consideră că vectorul x este fixat, existând o realizare a datelor, pentru care se efectuează testul Pentru demonstrație aplicăm metoda multiplicatorilor lui Lagrange, pentru a maximiza probabilitatea de detecție, cu constrângerea impusă de valoarea admisă pentru probabilitatea alarmei false Se construiește lagrangeanul F ⎡⎤ =+ −= + − αλαxx xxHH ()()() ⎢⎥ DFAFP P p d p dλ ; 10; ∫∫ ⎥ 11⎢⎣⎦ sau λα=+ −xxxHH ⎡⎤ ()() Fp p dλ ; 10; ∫ ⎣⎦ 1 Pentru a maximiza F, vom include în regiunea critică doar acei vectori x pentru care integrandul este pozitiv; integranzii pozitivi cresc valoarea integralei, în timp ce integranzii negativi scad valoarea integralei Prin urmare: 13 ∈+>xxxHH ()() 110; ppλ ; 0 Deoarece densitățile de probabilitate sunt pozitive inegalitatea de mai sus se pune sub forma () pxH 1; ∈>−x 1 λ () pxH 0; =-λ>0 sau, cu notația γ () pxH 1; ∈>x 1 γ () pxH 0; Pragul testului satisface ecuația de constrângere deja amintită γ===>xx xxH ()() {} FAPp d I Lα 0; ; unde : ∫ I *** Vom relua exemplul introductiv Raportul de plauzibilitate cerut de teorema NP se compară cu pragul testului, γși 2 − 11x⎧⎫ () ⎪⎪ exp− ⎨⎬ pπ 22; xH () ⎩⎭ 1⎪⎪ xH : Lγ ; == > () 12decidemdacă xH () () 0; 1px⎧⎫ ⎪⎪ exp− ⎨⎬ 22π ⎪⎪ ⎩⎭ După simpificări, condiția în care decidem că ipoteza H1 este adevărată, devine exp 1 2xγ−> {} sau, după logaritmare (care este posibilă, deoarece γ>0)și rearanjare ′ ln 1 2xγγ>+= Pragul testului se determină din ecuația de constrângere =−>==H {} {} 0 001FAPP xγα0exp 1 2 ; sau din echivalenta ei ′ =>==H {} 0 001FAPPxγα0 ; Repartiția eșantionului x în ipoteza nulă este 2 ⎫ 1x⎧ () ⎪⎪ H () ⎨⎬ exppx=−0 ; 22π ⎪⎪ ⎩⎭ așa că ecuația de constrângere a valorii probabilității alarmei false devine ∞ 21 ⎡⎤ =− == exp 2 ( ) 0 001PxdxQγ () FA′ ∫ ⎣⎦ 2π γ ′ ’=3 09 Probabilitatea detecției corecte este cea Se determină valoarea pragului,γ anterior determinată, de 0 0183 14 În cele ce urmează vom utiliza frecvent relația ′ μ− ⎛⎞ 2γ ′ >= dacă: xPxQγμσ∼N, , atunci: {} () ⎜⎟ σ ⎝⎠ Demonstrarea ei este un simplu exercițiu Expresia densității de probabilitate a unei variabile aleatoare normale este 2 ⎧⎫ − () 1xμ ⎪⎪ () ⎨⎬ 2exppx=− 22σ πσ ⎪⎪ ⎩⎭ ’ cu relația Se calculează probabilitatea ca variabila normală să depășească pragul γ ()/2uxμσ=− 2∞∞ ⎧⎫ − ⎧⎫ ′ () − 11xuμγμ ⎪⎪ ⎛⎞ ′ >= − − = = {} ⎨⎬ ⎨⎬ ⎜⎟ 2exp expPx dx du Qγ ∫∫ σσ πσ π ⎝⎠ ⎩⎭ − μσ ()/2222γγ ⎪⎪ ′′ ⎩⎭ exemplu în care analiză***Vom aborda un m detecția unei componente continue, din zgomot alb, gaussian Definim testul de ipoteze ipoteze sub forma =H 0: [ ] [ ] 0,1, , -1xn wn n N= += > H 1: [ ] [ ] 0,1, , -1 si 0xn A wn n N A= 22 CI∼ HH N () 01 si : [ ] 0, wuwnσσ = Vectorul de date este compus din cele N eșantioane de semnal măsurabil T [ 1]xx xN=−x [] Cele două ipoteze sunt caracterizate de medii diferite ale datelor T μ0H [] 0: 0 0 0== T μ1H [] 1: AAA A== Matricea de covarianță a zgomotului și deci și a datelor, în acest exemplu, este aceeași, în ambele ipoteze Rezultă că repartițiile vectorilor de date, în cele două ipoteze sunt 2 x0I∼HN () 0: , uσ 2 N x1I∼H () 1: , uAσ Testul de ipoteze este echivalent cu testul de parametri μ0H 0: = μ1H 1: A= rată (semnalul util e prezent) dacă Detectorul va decide că H1 este ipoteza adevă adică raportul de plauzibilitate depășește pragul γ, 1N− 211⎧⎫ () ⎨⎬ ∑ /22exp [ ]Nxn A−− nσ0222; = ⎩⎭ ()xH () 1pπσ x == > () 1NLγ xH 2;11p− () 0⎧⎫ () ⎨⎬ ∑ /22exp [ ]Nxn− = 0222nσ ⎩⎭ πσ () După simpilficare rezultă 11NN−− ⎡⎤ 21⎧⎫ −−+> ⎨⎬ ∑∑ 2exp 2 [ ]Axn Aγ ⎢⎥ 002nnσ== ⎣⎦ ⎩⎭ Se logaritmează 1N− 21⎡⎤ −− +> ∑ 22[] lnAxnNAγ ⎢⎥ 02nσ= ⎣⎦ 15 și se obține 12N− ANA NA >+ > ∑ 22 2[ ] ln se împarte cu 0xnγ σ 02nσσ= Obținem, în final, 12N− 1Aσ ′ γ=>+= [] lnxxnγ ∑ 02nNNA= În acest exemplu, detectorul NP compară media eșantion cu un prag Se știe deja că media eșantion este o estimare a componentei continue, A în cazul de față Dacă media eșantion are o valoare pozitivă, suficient de mare, probabil că este prezentă A>0 ’ se controlează probabilitățile alarmei false și de detecție corectă Știm Ajustând pragul γ că statistica Tx=x () are o repartiție gaussiană Drept urmare, vom determina mediile și dispersiile ei, în cele două ipoteze Avem 11NN−− 11⎧⎫ xH = (){} {} ⎨⎬ ∑∑ [] [] 0ET E wn Ewn==0; 00nnNN== ⎩⎭ 112NN−− 11σ ⎧⎫ xH = {} () {} ⎨⎬ ∑∑ [] []Disp T Disp w n Disp w n==02; 00nnNN N== ⎩⎭ 11NN−− 11⎧⎫ xH =+= ()(){} {} ⎨⎬ ∑∑ [] []ET E Awn EAwn A=+1; ⎩⎭ 00nnNN== 12N− 1σ ⎧⎫ xxHH = () () {}{} ⎨⎬ ∑ ; 10; [ ]Disp T Disp T Disp w n== 0nNN= ⎩⎭ În consecință, reprtiția statisticii T, în cele două ipoteze, este 2 ⎧ HN () Nσ în ipoteza 00, ⎪ T⎨ x∼ () 2 N H () ⎪ ANσ în ipoteza 1, ⎩ Aplicând relația roșie putem determina probabilitățile alarmei false și detecției corecte ⎛⎞ ′ 0γ − ′ α =>= = xH ⎜⎟ () {} FAPPT Qγ 0; 2⎜⎟ Nσ ⎝⎠ ⎛⎞ ′ Aγ − ′ xH =>= ⎜⎟ () {} DPPT Qγ 1; 2⎜⎟ Nσ ⎝⎠ Din prima relație se determină valoarea pragului, cu expresia 2 − 1σ ′ = () FAQPγ N Dacă substituim expresia pragului în relația de calcul a probabilității de detecție rezultă 22 ⎛⎞ 11NANA−− =− = ⎜⎟ ()() () DFA FAPQQP QQP ENRENR=−22; ⎜⎟ σσ ⎝⎠ Cu ENR am notat raportul dintre energia semnalului și puterea zgomotului (energy to noise ratio) În literatura de specialitate, de cele mai multe ori, acest raport este numit SNR 16 În figură sunt trasate curbele dependenței probabilității de detecție corectă, ca funcție de ENR, exprimat în dB, având ca și parametru de curbă probabilitatea alarmei false 1− () () DFAPQQ P ENR=− Admițînd α=0 001,la o valoare ENR de 14 dB (un raport de aproximativ 25 12) se obține o probabilitate de detecție corectă de 0 975, adică 97 5% Pentru ENR=20dB (raport de 100) 10− 10 0 9999PPα== ⇒ = FA D Exemplul analizat este o ilustrare a unei probleme mai generale de testare a ipotezelor, numită și problema repartițiilor gaussiene cu valori medii diferite În esență, se calculează valoarea unei statistici, T, și se compară cu un prag, γ ’ Dacă T depășește pragul, se declară adevărată ipoteza H1, altfel, H0 În ambele ipoteze avem 2 ⎧ HN () μσ00, în ipoteza ⎪ > μ∼ Tμ 10; 2⎨ H N () ⎪ μσ11, în ipoteza ⎩ În exemplul discutat, statistica T a fost media eșantion În acest caz, mai general, probabilitatea alarmei false este ′ − ⎛⎞ 0γμ ′ α H => = = {} FAPPT Qγ 0; ⎜⎟ σ ⎝⎠ ecuație din care rezultă valoarea pragului testului în această abordare ′ μ − 110γ−− ′ σ =⇒=+ () () FA FAQP QPγμ 0 σ Expresia valorii pragului testului se substituie în relația probabilității de detecție 0> ⎛⎞ ′ ⎜⎟ μμμ −− ⎛⎞ 1110γ− ′ H => = = − {}() DFAPPT Q QQPγ 1; ⎜⎟ ⎜⎟ σσ ⎝⎠ ⎜⎟ ⎝⎠ 17 Definim coeficientul de deflexie cu relația 2 − () 210μμ 2d= σ care înlocuită în expresia probabilității de detecție dă 12− () ) ( DFAPQQP d=− Q(x) fiind o funcție descrescătoare, probabilitatea de detecție crește, odată cu creșterea coeficientului de deflexie 18